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ABSTRACT 

We present a sample of 75 extinction curves derived from FUSE far-ultraviolet spectra supplemented 
by existing IUE spectra. The extinction curves were created using the standard pair method based 
on a new set of dereddened FUSE+IUE comparison stars. Molecular hydrogen absorption features 
were removed using individualized Hi models for each sightline. The general shape of the FUSE 
extinction (8.4 /im" 1 < A -1 < 11 /1m" 1 ) was found to be broadly consistent with extrapolations 
from the IUE extinction (3.3 /im -1 < A -1 < 8.6 /im" 1 ) curve. Significant differences were seen in 
the strength of the far-UV rise and the width of the 2175 A bump. All the FUSE+IUE extinction 
curves had positive far-UV slopes giving no indication that the far-UV rise was turning over at the 
shortest wavelengths. The dependence of A(X)/A(V) versus i?(U) _1 in the far-UV using the sightlines 
in our sample was found to be stronger than tentatively indicated by previous work. We present an 
updated R(V) dependent relationship for the full UV wavelength range (3.3 /im -1 < A -1 < 11 ^ixnT 1 ). 
Finally, we searched for discrete absorption features in the far-ultraviolet. We found a 3cr upper limit 
of ^Q.12A(V) on features with a resolution of 250 (~4 A width) and 3<r upper limits of ~0.15A(V) 
for A -1 < 9.6 iiraT 1 and ^0.68A(V) for A -1 > 9.6 /im -1 on features with a resolution of 10 4 (~0.1 A 
width). 

Subject headings: dust, extinction 
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1. INTRODUCTION 

Extinction curves are one of the cornerstones of our 
understanding of dust in interstellar space. Extinction 
curves probe the sum of the absorption and scattering of 
dust and their determination is observationally straight- 
forward. Most extincti on curves are measured using the 
standard pair method (|Stecherlll965l ; IMassa et al.lll983h 
where the measurements of a reddened star are ratioed 
to those of an unreddened star of the same spectral 
type. The importance of measuring extinction curves 
over the widest wavelength range possible is attested to 
by their use as a fundamental constraint in models of 
dust grains (IWeingartner fc Draine1l2001t IClavton et al.1 
l2003bl : IZubko et al.ll2004D . Measuring extinction curves 
is also important for purely empirical reasons in order to 
allow for the proper accounting for the effects of inter- 
stellar dust in the study of astrophysical objects. 

Structure in interstellar extinction curves gives di- 
rect evidence of the dust grain compositions. In 
the ultraviolet, the one very ob vious discrete featur e 
is the 2175 A extinction bump (|Stecherl I1965L I1969D . 
The 2175 A bump has been observed to have nearly 
constant ce ntral wavelength, a width varying from 
280-660 A (iFitzpatrick fc Massal Il986t IValencic et~ai1 
120041: IFitzpatrick fe Massal 120071). and has been at- 
tributed to small graphit e grains (jStecher fc Dondll965l : 
iDraine fc~M alhotra 1993]). There have been no other dis- 
crete features found in the ultrav iolet with fairly sensi- 
tive limits placed from searches (jClayton et al.l [2003a) 
for counterparts of the optical/ near- infrared Diffuse In- 
terstellar Bands (;DIBs. lMerrTi|[il3l:lHerbr3 ri995)^ The 



other main structure seen in the ultraviolet extinction 
curve is the far-ultraviolet (far-UV) rise. This rise is 
characte rized by a constant shape between 1700 and 
1175 A (|Fitzpatrick fc Massalll988l ) and has been mod- 
eled as the wing o f a feature like the 2175 A bump, bu t 
centered at 715 A (jJoblin et al.lll992tlLi"fc Drainell2001h . 

The structure in the ultraviolet between 1175 and 
3300 A has been well studied using spectra taken 
with the International Ultraviolet Explore r (IUE) 
(IFitzpatrick fc Massalfl990l : IValencic et al.ll2004h and the 
Space Tele scope Spectro g raph o n the Hubble Space 
Telescope (jClavton et al.l l2003al ). The structure in 
this wavelength range is well approximated by the 
IFitzpatrick fc Massal |l990f l (hereafter FM90) parame- 
terization with the curve being the combination of a lin- 
ear term, a Drude profile, and a cubic describing the far- 
UV rise. In addition, th e average beha v ior of ultraviolet 
extinction was found bv lCardelli et al.l (|1989f ) (hereafter 
CCM89) to correlate with R{V) = A(V)/E(B - V) (a 
rough measure of the average dust grain size). More re- 
cently, IFitzpatrick fc Massal (|2007l ) have proposed a re- 
vised parameterization that includes an additional pa- 
rameter and simplifies the far-UV curvature term to just 
a quadratic. 

The dust extinction in the far-UV region between 912 
and 1190 A has onl y been studied fo r a handful of Milky 
Way sightlines (see lSofia et al.l l2005. for a review). Our 
understanding of the far-UV extinction and the nature of 
the dust materials responsible for it is, necessarily, tenta- 
tive. The many spectra of hot stars taken w ith the Far- 
Ultraviolet Spectroscopic Explorer (FUSE, iMoos et al.l 
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TABLE 1 
Comparison Stars 



name 


SpType 


V 


A(V) 


R(V) 


Ci 


C 2 


c 3 


C 4 


x 


7 


V, main sequence 


HD091824 


06V 


8.14 


0.60 ± 0.06 


3.10 ± 0.31 


-0.81 ± 0.41 


0.90 ± 0.22 


3.69 ± 1.11 


0.43 ± 0.06 


4.62 ± 0.09 


1.00 ± 0.10 


HD093028 


08V 


8.37 


0.80 ± 0.12 


3.20 ± 0.48 


0.38 ± 0.19 


0.61 ± 0.15 


2.07 ± 0.41 


0.29 ± 0.06 


4.62 ± 0.18 


0.87 ± 0.09 


HD097471 


BOV 


9.30 


0.90 ± 0.09 


3.10 ± 0.31 


-0.93 ± 0.47 


0.88 ± 0.18 


1.83 ± 0.55 


0.28 ± 0.06 


4.62 ± 0.09 


0.85 ± 0.09 


BD+52°3210 


B1V 


10.69 


0.85 ± 0.09 


3.50 ± 0.35 


-1.67 ± 0.42 


0.99 ± 0.15 


3.43 ±0.69 


0.25 ± 0.06 


4.60 ± 0.09 


0.91 ± 0.09 


BD+32°270 


B2V 


10.29 


0.10 ± 0.04 


3.10 ± 1.24 


-0.24 ± 0.12 


0.74 ±0.37 


2.00 ± 1.00 


0.82 ± 0.41 


4.59 ± 0.46 


1.00 ± 0.20 


HD051013 


B3V 


8.81 


0.10 ± 0.04 


3.10 ± 1.55 


-0.24 ± 0.12 


0.74 ±0.37 


2.00 ± 1.00 


0.42 ± 0.21 


4.59 ± 0.46 


1.00 ± 0.20 


HD037332 


B4V 


7.62 


0.10 ± 0.02 


3.10 ± 0.60 


-0.24 ± 0.12 


0.94 ± 0.47 


3.56 ± 1.78 


2.00 ± 0.50 


4.59 ± 0.46 


1.00 ± 0.20 


HD037525 


B5V 


8.08 


0.18 ± 0.09 


3.10 ± 1.24 


-0.24 ± 0.12 


0.74 ± 0.37 


3.56 ± 1.78 


0.52 ± 0.26 


4.59 ± 0.18 


1.00 ± 0.20 


III, giant 


HD116852 


09III 


8.47 


0.63 ± 0.06 


3.01 ± 0.60 


-0.65 ± 0.33 


0.70 ± 0.10 


1.41 ± 0.28 


0.28 ± 0.07 


4.56 ± 0.18 


0.72 ± 0.07 


HD104705 


BOIII 


7.76 


1.06 ± 0.11 


3.31 ± 0.33 


0.27 ± 0.13 


0.49 ± 0.10 


2.01 ± 0.40 


0.30 ± 0.03 


4.56 ± 0.09 


0.82 ± 0.08 


HD 172 140 


B0.5III 


9.94 


0.56 ± 0.11 


2.34 ± 0.47 


-0.37 ± 0.18 


0.61 ± 0.15 


2.58 ± 0.52 


0.37 ±0.06 


4.59 ± 0.09 


0.89 ± 0.09 


HD 114444 


B2III 


10.31 


0.51 ± 0.10 


2.65 ± 0.53 


-0.65 ±0.32 


0.40 ± 0.10 


2.56 ± 0.51 


0.31 ± 0.08 


4.61 ± 0.09 


0.82 ± 0.08 


HD235874 


B3III 


9.64 


0.58 ± 0.06 


2.93 ± 0.29 


-2.63 ± 0.66 


1.11 ± 0.11 


3.26 ± 0.65 


0.26 ± 0.06 


4.68 ± 0.09 


0.97 ± 0.10 



I, supcrgiant 



HD210809 
HD091983 
HD094493 
HD100276 
HD075309 



091b 
09.51b 
B0.5Iab 
Bllb 
B2Ib 



7.54 
8.58 
7.27 
7.22 
7.85 



0.98 ± 0.05 
1.10 ± 0.11 
0.70 ± 0.04 
0.83 ± 0.08 
0.64 ± 0.06 



3.10 ± 0.31 
3.10 ± 0.23 
3.30 ± 0.33 
3.50 ± 0.35 
3.53 ± 0.71 



-1.23 ± 0.31 
0.35 ± 0.17 
-1.02 ± 0.51 
-1.48 ±0.37 
0.39 ± 0.20 



0.98 
0.61 
0.95 
0.92 
0.53 



± 0.05 
± 0.09 
± 0.19 
± 0.09 
± 0.13 



2.58 ± 0.26 
2.50 ± 0.50 
3.14 ± 0.63 
3.30 ± 0.66 

2.59 ± 0.52 



0.45 ± 0.02 
0.40 ± 0.04 
0.56 ± 0.08 
0.26 ± 0.06 
0.30 ± 0.08 



4.60 ± 0.09 

4.61 ± 0.09 
4.60 ± 0.09 

4.59 ± 0.09 

4.60 ± 0.09 



1.00 ± 0.10 
0.98 ± 0.10 
0.93 ± 0.09 
0.95 ± 0.10 
0.76 ± 0.08 



2000) allow for an extensive study of the extinction curve 
in the far-UV region and motivate this analysis. The aim 
is to refine our understanding of the far-UV extinction 
using a statistically significant sample of FUSE extinc- 
tion curves. Specific ally, we will investigate how well the 
FM90, CCM89, and lFitzpatrick fc Massal (|2007t) param- 
eterizations extrapolate to the far-UV and refine them if 
necessary. Also, we will see if the far-UV region contains 
any significant structure in addition to the well-known 
far-UV rise. 

This study is the third paper in a series of papers in- 
vestigating FUSE fa r-UV extinction curves. The first 
paper in this series (|Sofia et al.l l2005f) presented a pre- 
liminary study of far-UV extinction curves along nine 
Milky Way sightlines. The second paper in this series 
(|Cartledge et al.ll2005f ) presented a study of nine far-UV 
extinction curves in the Magellanic Clouds. This third 
paper presents far-UV extinction curves for 75 sightlines 
and benefits from a more mature FUSE calibration, a 
better understanding of correcting for the H2 absorption, 
and a more complete set of FUSE comparison stars . The 
fourth paper in this series ICartledge et al.l (|2009f ) will 
study the correlations between the gas properties (e.g., 
H2, & H I) and dust extinction. 



corrections for the copious H2 absorption seen through- 
out the FUSE spectral region. This was done by deter- 
mining the physical parameters of the H2 absorbing gas 
and generating a model spectrum of the H2 absorption. 
The details of t he H2 and H I modeling are given in 
ICartledge et al.l (|2009h . When correcting a FUSE spec- 
trum for H2 and H I absorption, any point which is ab- 
sorbed by more than 70% was excluded from use in gen- 
erating an extinction curve, to minimize the effects of 
the uncertainties in the H2 and H I corrections on the 
resulting extinction curves. The residual effects of the 
H2 and H I absorption on the final extinction curves are 
discussed in 33.41 

The complimentary IUE spectra (1150 - 3225 A) were 
taken from the IUE archive following VCG04. Optical 
and near-infrared photometry was accumulated from the 
literature also following VCG04. The spectral overlap 
between the IUE and FUSE observations was used to 
derive a multiplicative correction to the FUSE spectra 
to put them on the same flux scale as the IUE spectra. 
The corrections are generally small with maximum cor- 
rections up to 25%. Finally, the FUSE and IUE spectra 
were rebinned to a common resolution of 250 which is 
sufficient for most extinction curve work. 



2. DATA 

The sample of stars for this paper was determined 
by searching the F USE archive for stars from the 
iValencic et atl (l200l (hereafter VGC04) sample of IUE 
extinction curves that had good quality FUSE spectra. 
IUE extinction curves are required because the large scale 
extinction structure in the FUSE wavelength range (far- 
UV rise) starts in the IUE s pectral range around 1700 A 
(|Fitzpatrick &: Massalfl9 88). As a side benefit, the exis- 
tence of IUE extinction curves for all the sightlines stud- 
ied in this paper ensures that these sightlines are suitable 
for extinction curve determinations. The final sample for 
this paper includes 75 sightlines. 

The FUSE observations (912 - 1190 A) were ex- 
tracted from the archive and reduced using CalFUSE 
v3.0. The obse r vation s used for each star are given in 
ICartledge etaLl (|2009l ). The next step was to generate 



2.1. New Comparison Stars 

Calculating an extinction curve for the sightline to- 
wards a reddened star requires an unreddened compar- 
ison star of the same spectral type. For IUE obser- 
vations, there exists a good set of bright, lightly red- 
dened comparison sta rs which have been carefully dered- 
dened (|Cardelli et al.lll992l) . Ideally, FUSE observations 
of these same IUE comparison stars would be used to ex- 
tend their spectra to far-UV wavelengths. Unfortunately, 
many of the IUE comparison stars are too bright to be 
observable with FUSE. A new set of fainter comparison 
stars that are observable with FUSE is required. 

We have generated the FUSE set of comparison stars 
by picking the least reddened star of each spectral type 
that has both FUSE and IUE observations. The new 
comparison star sample is listed in Table[T]along with the 
parameters used to deredden their IUE and FUSE spec- 



FUSE Extinction Curves 
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Fig. 1. — The FUSE and IUE spectra of the comparison stars are plotted. All the spectra are plotted at a common spectral resolution 
of 250. The plots on the left give the observed spectra and the same spectra which have been dereddened and corrected for H I and H2 
absorption are on the right. The spectra are normalized to one between 4 and 5 (im" 1 and offset from each other by 0.5 in log space. 
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tra. The H2 and H I model parameters used to remove 
the a ffects of H2 absorption are given by (jCartledge et al.1 
2009). Given the FUSE brightness limit, these stars 
are necessarily more distant and reddened than the IUE 
comparison star sample. We have carefully dereddened 
these new comparis on stars following the methods of 
ICardelli et al.l (|1992f ). The deredd ening was done start- 
ing w ith E{B - V), A(V), and IFitzpatrick fc Massal 
( 1990) parameters taken from VCG04. Thus, these new 
comparison stars are directly bootstrapped off of the ex- 
isting IUE comparison stars. The use of FM parameters 
determined in the IUE spectral range to deredden FUSE 
spectra is supported by earlier work on FUSE extinc- 
tion curves showing that the extrapolation of the FM fits 
to th e FUSE wavelength range is reasonable l) Sofia et al.l 
I2005T ) . We have checked that this extrapol ation is not af- 
fecting our results in H3.3I As was done bv lCardelli et al.1 
(1992), the dereddening parameters were then manually 
tweaked to produce dereddened spectra that lacked ob- 
vious dust extinction and had a smooth progression be- 
tween spectral types. This is illustrated in Fig. [T] where 
the observed and dereddened spectra for each class of 
standards (main sequence, giants, and supergiants) are 
shown. The uncertainties on the dereddening parameters 
was estimated by changing each parameter until it was 
clearly incorrect. The resulting uncertainties are given in 
Table [1] and are assumed to be one sigma uncertainties 
to be conservative. 

An alternative to using observed comparison stars is to 
use stellar atmospher e models. This is an approach that 
has been explored by IFitzpatrick fc Massal (|1999L l2005f ) 
and found to produce similar results to using observed 
comparison stars. The strength of using model atmo- 
spheres is that it removes the need to deredden observed 
spectra and should provide closer spectral matches. The 
weakness is that it relies on the accuracy of the stel- 
lar atmosphere models and the absolute calibration of 
the spectra. Observed comparison stars only rely on the 
relative calibration of the spectra and, by definition, in- 
clude the same physics as the reddened stars. The use 
of observed comparison stars does inject some additional 
uncertainty due to the lack of "perfect" matches to the 
reddened star spectral types. In the end, using observed 
or model comparison stars represent complimentary ap- 
proaches, each with strengths and weaknesses. The fact 
that they produce similar results provides strong evi- 
dence that the resulting extinction curves derived from 
either method are correct. 



2.2. Extinction Curve Calculation 

The extinction curves were calculated using the stan- 
dard pair method (|Stecherlll96l iMassa et al.lll983[ ) . Ba- 
sically, the ratio of the fluxes of the reddened and com- 
parison stars (both with the same spectral type) gives a 
direct measurement of the dust extinction towards the 
reddened star. Ideally, the distances to the reddened 
and comparison stars would be known to high accuracy, 
allowing for an absolute measurement of the dust extinc- 
tion. Unfortunately, the distances to the hot, early-type 
stars used for ultraviolet dust extinction curve measure- 
ments are rarely, if ever, known to high enough accuracy. 
It is possible then to make a differential measurement of 
the difference between the colors of reddened and com- 



parison stars. The basic measurement is 

E(X -V) = m(A - V) r - m(A - V) c (1) 

where r and c subscripts refer to the reddened and com- 
parison stars, respectively. This measurement is the un- 
normalized extinction curve. The use of colors w.r.t. the 
V band removes the unknown or uncertain distance to 
each of the stars from the measurement. In order to 
make comparisons between extinction curves measured 
along different sightlincs, the curve needs to be normal- 
ized to a measurement of the total column of dust along 
the sightlinc. The most often used normalization is to di- 
vide E(X — V) by E(B — V). Yet a more direct measure 
of the dust properties is A(X)/A(V). 

There are two avenues to determining A(X) /A(V) from 
the basic E(X — V) measurement. The usual method is to 
determine the conversion factor R(V) [— A(V)/E(B — 
V)} by extrapolating the E(X - V) /E(B - V) curve at 
JHK wavelengths to infinite wavelength using an as- 
sumed form of the A(X)/A(V) curve at wavelengths 
> 1 /im. The conversion is then 



1 



A{\) 

A(V) R(V) 



E(X - V) 



E(B - V) 



1 



(2) 



A more direct way is to extrapolate the E(X — V) curve 
to infinite wavelength using the same assumptions [JHK 
wavelengths and an assumed A(X) / A(V) curve] to deter- 
mine A(V). The conversion is then 



A(X) _ E(X - V) 
W)~ A{V) 



1. 



(3) 



Functionally, both methods of determining A(X)/A{V) 
are equivalent. 

The use of A(X)/A(V) extinction curves is preferred 
over performing our analysis on E(X — V)/E(B — V) 
extinction curves as A(X)/A(V) is the more fundamen- 
tal measurement of the properties of dust grains. In 
addition, the A(X)/A{V) curve is less affected by sys- 
tematic uncertainties introduced by the normalization. 
This is not surprising as the fractional uncertainty on 
A(V) is lower (on average 3x) than the fractional un- 
certainty on E(B — V). This is due to the quality of 
the 2MASS observations that provide the JHK photom- 
etry and the fact that the A(V) measurement is based 
on (effectively) the average of the JHK measurements of 
the extinction curve. Of course, the A(X)/A(V) curve 
does include the assumption that the > 1 /jm curve can 
be extrapolated accurately to derive A(V). This ap- 
pears to be a reasonable assumption for the diffuse ISM 
(|Martin fc Whitte~tlll990l ). but may not be for sightlines 
that probe dense IS M regions (|Indebetouw et all 120051 : 
iFlahertv et~aTll2007fh Most (if not all) of the sightlines 
in this pap e r qua lify as diffuse sightlines. We use the 
iRieke et al.1 (|198 9) infrared extinction curve to extrapo- 
late the JHK E(X — V) curves to infinite w avelength to 
derive A(V) values. ThelRieke et al.l (119891) work refines 
the results of lRieke fc Lebofskvl 1)19851 ) . 

One useful benefit of directly determining A(X)/A(V) 
from E(X — V) is that this simplifies the calculation 
of the uncertainty on the normalized extinction curve. 
Previously, we would calculate E(X — V)/E(B — V) 
curve, determin e R(V), and then renorm alize the curve 
to A(X)/A(V) dGordon fc Clavtonl [l998h . Propagating 
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FlG. 2. — The extinction curves are plotted for the entire sample 
by measured R(V) values and offset by a constant value from each 
the right of each curve. 
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along with the FM90 fit to each curve. The curves are plotted sorted 
other. The reddened star name and measured R(V) value is given to 
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TABLE 2 
Extinction Curve Details 



name SpTypc comparison V A(V) a E(B-V) a R(V) 



BD+35°4258 


B0.5V 


HD097471 


9 


.41 





.95 


± 





.04 ± 





.03 





.31 


± 





.053 


± 


.037 


3 


.06 


± 


.380 


± 


Oil 


BD+53°820 


BOIV 


HD104705 


9 


.95 


1. 


.20 


± 





.04 ± 





03 





.40 


± 





.047 


±0 


.035 


3 


.01 


±0 


.249 


±0 


.012 


BD+56°524 


B1V 


BD+52°3210 


9 


.75 


1. 


.80 


± 





.04 ± 





.02 





.60 


± 





.044 


± 


.023 


2 


.98 


± 


.191 


± 


.009 


HD001383 


Bill 


HD100276 


7 


.63 


1. 


.36 


± 





.04 ± 





.02 





.45 


± 





.035 


± 


.029 


3 


.04 


± 
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a The quantititics are given as value ± random uncertainty ± systematic uncertainty. 
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Fig. 3. — The A(V) and R(V) values for each extinction curve are 
plotted. The histograms of both quantities also help characterize 
the sample. 

uncertainties using these steps adds the uncertainty in 
E(B — V) twice in the final A(X)/A(V) curve when, in 
fact, the E(B — V) uncertainty does not contribute to the 
final uncertainties in A(X)/A(V) at all. This resulted in 
the unc ertainties we have quoted in some of our p revious 
papers (jGordon et alj l2003t IValencic et all I2004D being 
overestimates. The correct A(X)/A(V) uncertainties are 



A(V) 



'_(A{X) 
\A{V) 

+ ^(A) 2 




o[E{\-V)\ X' 
E(X-V) J 



A(V) 



where 



a[E(X-V)] = \/<r[m(X)r} 2 + a[V r } 2 + a[m(X) c } 2 + a[V c } 2 

(5) 

and CTd(A) is the uncertainty due to dereddening the com- 
parison stars f ^2.ip . The value of aa(X) for each curve 
was determined using a Monte Carlo method where the 
dereddening is done using parameters picked from the 
distribution allowed by the uncertainties on the dered- 
dening parameters and computing the resulting uncer- 
tainty on each A{X)/A(V) point. As can be seen from 
these equations, the A(X)/A(V) uncertainties can be di- 
vided into random (varying from point to point) and sys- 
tematic (causing the entire curve to move up or down) 
terms. The random uncertainties are due to the <r[m(A) r ] 
and cr[m(A) r ] terms. The systematic uncertainties are 
due to the a[V r ], a[V c ], a[A(V)], and <Jd{X) terms. 

The FUSE+IUE extinction curves for all the stars in 
our sample are shown in Fig. [2] The comparison stars 
used along with details of each extinction curve are given 
in TableO The quantities are given as value ± random ± 
systematic uncertainties. The comparison star used for 
each reddened star was picked to be the closest in the 2- 
dimensional space defined by the MK spectral type. We 
rejected any extinction curves that had a reddened star 
to comparison spectral-type mismatch larger than one 




hd062542; R(V)=3.1' 



1/> [Mm" 1 ] 

FlG. 4. — The updated ver sions of 6 of the extinction curves 
presented in[Spfia ct al. (2005) arc shown. 

2 \ 

temperature or luminosity subclass to account for spec- 
/ tral type uncertainties. The A(V) and R(V) distribution 
(4)°f our sample are shown in Fig. [3l 

3. RESULTS 
3.1. Comparison to \Sofia et al\ 11200 A ) 

The preliminary work o n this topic u sing F USE obser- 
vations was presented in ISofia et all (|2005f ). The new 
extinction curves presented in this paper improve on 
this previous work in a number of ways. First, the 
calibration pipeline version used was v3.0 (compared 
to v2.2.1) which significantly improved the reduction 
of faint sources, especially at the shorter FUSE wave- 
lengths. Second, only sightlines where the H2 modeling 
included populations beyond J=l were included. Third, 
the set of comparison stars included more spectral and 
luminosity types and the comparison star spectra were 
carefully dereddened. Finally, our new work imposed 
more stringent data quality cuts and a requirement for 
close spectral type matches between reddened and com- 
parison stars. T hese changes resulted in three of the 
I Sofia et ail (|2005f ) sightlines being rejected from our sam- 
ple. These three stars were HD 167971 (poor H2 model), 
HD 210121 (no good comparison star), and HD 239682 
(poor FUSE and IU E data) 

The six stars from lSofia et al.l (|2005| ) remaining in our 
sample are shown in Fig. [4] These e xtinction cu r ves ar e 
similar to those shown in Fig. 2 of So fia et ail ([2005) . 
They extend to shorter wavelengths and are generally 
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Fig. 5.— The FM90 coefficients fit to the IUE+FUSE data are 
without the FUSE data. In determining the fitting uncertainties for 
systematic uncertainties affect the IUE+FUSE and IUE only curves 

less noisy. 



3.2. FM Fits 

Each A(X)/A(V) curve was fit with the FM90 
arameterization of the UV extinction curve 
Fitzpatrick fc Massal 1 19861 I1988L Il990t ). This pa- 
rameterization was developed to describe the extinction 
curve in the IUE spectral range (1150 - 3300 A). The 
FM90 function, modified to fit A(X) /A(V) curves, is: 

(6) 



plotted versus the difference between these coefficients fit with and 
this comparison, only the random uncertainties were considered. The 
in the same manner. 

where x = l/\\fj,m~ 1 ], the 2175 A bump is represented 
by the Drude term 

2 

° {X > 7 ' Xo) = (a; 2 - xlY + (x 7 ) 2 ' (?) 

and the far-UV curvature (for x > 5.9 /mi -1 ) is given by 

F{x) = 0.5392(a; - 5.9) 2 + 0.05644(x - 5.9) 3 (8) 

These 6 parameters are not uncorrelated and deter- 
mining the best fit requires iteration to converge to a 
stable solution. We have found that a fitting algo- 
rithm that alternates between fitting the four coefficients 



A(V) (jMV) n A(V) 



A(V). 



, 7, and x produces the 



FUSE Extinction Curves 
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.99 


± 





019 


± 





.012 


HD239729 


f) 


.87 


± 





049 


± 


0. 


.021 





.26 


± 


0. 


009 


± 





.008 


1 


16 


± 





.060 


± 





.029 


0. 


.18 


± 





.008 


± 





004 


4. 


60 


± 


0. 


.022 


± 





.009 


1 


.11 


± 


0. 


021 


± 





.011 


HD326329 


1 


.07 


± 





.074 


± 


0. 


.077 





.27 


± 


0. 


013 


± 





.020 


1. 


.20 


± 





.089 


± 





.093 


0. 


12 


± 





.009 


± 


0. 


009 


4. 


61 


± 


0. 


033 


± 


0. 


.016 





.95 


± 





024 


± 





.018 


HD332407 





.64 


± 


0. 


050 


± 





.025 





.31 


± 


(J 


.009 


± 





.011 


1. 


27 


it 





.067 


± 





.051 


0. 


10 


± 





.007 


± 





.004 


4. 


57 


± 





.024 


± 


0. 


010 


1 


.01 


± 


0. 


019 


± 





.011 



The quantitities are given as value it random uncertainty it systematic uncertainty. 
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most stable and best quality fits. As part of this fitting 
algorithm, we also have found that iterative sigma clip- 
ping is important to ensure that a small number of highly 
deviant points does not bias the fits. The relationship be- 
tween our version of the FM90 parameters and tho se fit 
to E{\ - V)/E(B - V) (|Fitzpatrick k Massalll99l is 

cf v) =d/R(V) + l (9) 

Cf y) =Q/i?(n for j = 2, 3, 4. (10) 

The uncertainties on the resulting FM90 parameters arc 
composed of random and systematic components. The 
random component is due to the uncertainties on the 
measurement of the ultraviolet fluxes and the system- 
atic component is due to uncertainties in the V band 
magnitudes and in the calculation of A(V). The random 
uncertainties were calculated using the Monte Carlo tech- 
nique generating trials near the best fit and determining 
when the resulting FM90 fit is within the errors using 
the F-test method. The random uncertainty is then one 
half the difference between the minimum and maximum 
fit values that satisfy the F-test criteria. The systematic 
uncertainties are determined by one half the difference of 
the FM90 fit coefficients determined from the two curves 
generated by adding and subtracting the systematic un- 
certainties from the curve. Finally, the random and sys- 
tematic uncertainties can be combined in quadrature to 
produce the final FM90 parameter uncertainties. 

The FM90 fit coefficients for the extinction curves are 
given in Table [3] with the value ± random ± systematic 
uncertainties. 

Fitting the IUE+FUSE extinction curve with the 
FM90 parameterization does produce different FM90 fit 
coefficients than fitting just the IUE data. The effect 
of adding the FUSE data on the FM90 fit coefficients is 
shown in Fig. [5] The FM90 coefficients that show signif- 
icant change are and 7. This is not a surprising 
result as the Cf^ V ' coefficient describes the strength of 
the far-UV rise and 7 is the width of the 2175 A bump. 
The FUSE data put much stronger constraints on the 
far-UV rise and this, in turn, affects the best fit width 
of the bump. The most significant result is that there is 

an average shift in the value of C± (V) of 0.0072 ± 0.0016 
which is a ~8% reduction in the average strength on the 
far-UV rise. None of the other parameters have signifi- 
cant shifts in their average values. Finally, the addition 
of the FUSE data removes the negative far-UV curva- 
ture strengths (C^ V ^) that were seen with fits to only 

the IUE data. 

Recently, iFitzpatrick k Massal (|2007f l presented a 
modification to the FM90 parameterization (hereafter 
called FM07). They introduced a 7th coefficient (C 5 ) to 
allow the wavelength where the far-UV curvature term is 
important to be a fitted parameter. They also modified 
the far-UV term to just be a quadratic (removing the cu- 
bic term). We have compared the fit x 2 values between 
FM90 and FM07 fits to the IUE+FUSE extinction curves 
in Fig. [6l On average, there is a preference for the FM07 
over the FM90 parameterization, but this is a weak pref- 
erence given that 40% of the sample is better fit with the 
FM90 parameterization. There are significant differences 
between some of the parameters (C3 , C^ V \ and 7) 
between the FM90 and FM07 fits, but it could be that 
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Fig. 6. — The fractional difference between the x 2 values for 
the FM90 and FM07 fits are shown. On average, the FM07 x' 2 is 
better, with the split of which fit produces the smaller x 2 being 46 
to 30 between FM07 and FM90, respectively. 

this is just a reflection of the non-orthogonal nature of 
the fit parameters and/or the result of adding an addi- 
tional free parameter. Since there is no strong reason to 
prefer the newer FM07 parameterization over the FM90 
parameterization, we use the FM90 parameterization for 
the majority of this paper as it has one fewer free param- 
eter and is easier to compare with previous work. 

3.3. R(V) Dependent Relationship 

The work of CCM89 demonstrated that on aver- 
age extinction curves follow a family of curves that 
can be described with one parameter. Real deviations 
from the CCM89 r e lation ship exist and were studied by 
IMathis k CardeiTil (fl992h and VCG04. CCM89 chose 
R(V) as the one parameter as it roughly measures the 
average grain size. The CCM89 relationship is an empir- 
ical relationship based on determining the linear corre- 
lation between A(X)/A(V) and i?(U) _1 at each wave- 
length with the resulting set of intercepts and slopes 
forming the the CCM89 R(V) dependent relationship. 
The CCM89 relationship w as constructed mainly using 
IFitzpatrick k Massal (|1988f ) fits to IUE extinction curves 
and broadband data in the near-infrared and optical of 
29 sightlines. This was supplemented with ANS UV pho- 
tometry and Copernicus far-UV spectra for a smaller 
number of sightlines. The CCM89 relationship in the 
far- UV (8-10 jum -1 ) is mainly based on extrapolating 
the IFitzpatrick k Massal (|1988f) fits and represents the 
the most uncertain portion of the CCM89 relationship. 
Overall, CCM89 is valid from 0.3-10 ^m" 1 (3.3-0.1 fan). 

Since CCM89, there have been two other efforts to 
refin e the R(V) dependent relationship in the ultravi- 
olet (lFitzpatricldll99fl IValencic et al.ll2004D . iFitzpatrickl 
(1999) (hereafter F99) used the FM90 set of IUE extinc- 
tion curves as the basis for their R(V) dependent rela- 
tionship. Their sample included 77 sightlines and the F99 
relationship in the ultraviolet is based on the correlation 
between C\ and C 2 FM90 coefficients and R{V)~ l (the 
other coefficients are held constant). The near-IR and 
optical regions arc fit using splines to account for the 
broad-band nature (and shifting effective wavelengths) 
of the data. The F99 relationship is valid from ^0.2-8.7 
^m- 1 (5.0-0.115 fim). 

VCG04 used a sample of 417 IUE extinction curves as 
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Fig. 7. — The relationship between A(\)/A(V) and R(V) 1 is shown for selected wavelengths using FUSE (left) and IUE (right). The 
relationships given by the CCM89 (A < 10 Atm -1 ), F99 (A < 8.7 /im" 1 ), and VCG04 (A < 8 /im" 1 ) fits are also shown, but only for 
wavelengths for which they are valid. The IUE plot is similar to that of Fig. la of CCM89. The F99 relationship is non-linear in this form 
as it is based on correlations between FM90 fit parameters and iJfV)" 1 . 




Fig. 8. — The intercept (a) and slope (b) coefficients to the fits of A(\) / A(V) versus R(V) 1 are plotted versus wavelength. Only fits 
using 40 or more points are plotted. The coefficients for the individual fits (A fits) and the CCM89-type functional fit to them are shown. 
The coefficients from CCM89 and VCG04 are also given for comparison. The regions around the wind lines (6.5 and 7.1 /im~ , open 
circles) are not used in the CCM89-like fit as they are almost never well matched by the comparison stars. 



the basis for their R(V) dependent relationship. They 
followed the methodology of CCM89 and fit A(X) /A(V) 
and i?(U) _1 at each wavelength. The resulting slopes 
and intercepts were fit with the same functional form as 
CCM89. This work concentrated on the IUE UV region 
and this is only valid between 3.3-8.0 /urn -1 (0.3-0.125 

Using our sample of 75 FUSE extinction curves we can 
investigate the R(V) dependent relationship in the far- 
UV (8.4-11 /im- 1 ). We start by following the CCM89 



methodology and perform fits of the data to 
A(x) , , b(x) 

where x = A . For consistency, we have determined 
the linear fit coefficients for both our FUSE and IUE ex- 
tinction curves as this allows us to directly compare our 
R(V) dependent relationship to the CCM89, F99, and 
VCG04. Unlike CCM89, we perform our linear fits on 
the actual extinction curves, not the FM90 fits. This 
directly connects our R(V) relationship to the data and 
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does not assume a functional form for the UV extinction 
curve. The linear fits are done using the 'fitexy' IDL pro- 
gram that takes into account uncertainties on both the 
independent and dependent variables. Another approach 
was taken in deriving the R{V) dependent relationship 
by F99 who correlated the FM90 parameters with R(V). 
We did a similar analysis and found, just as F99 did, that 
most of the dependence of the extinction curves on R(V) 
is due to correlations with C\ and C 2 . We choose to use 
the CCM89 methodology for our R(V) relationship as it 
imposes fewer assumptions about the detailed structure 
of the relationship. Both approaches (CCM89 or F99) 
produce very similar R(V) relationships. 

The fit relationships, fits, and comparisons to previ- 
ous R(V) dependent relationships are shown in Fig. [7l 
The relationship shown by our sample of sightlines is 
similar to those seen in previous works, but our fits dis- 
play a stronger dependence on i?(V) _1 . In addition, our 
R(V) dependent relationship is well defined all the way 
to 11 /im -1 . As can be seen from this figure, our R(V) 
dependent relationship is valid for R(V) values between 
2.5 and ~5. The intercept (a) and slope (b) values are 
shown as a function of wavelength in Fig. [8] along with 
the values from previous works. We find quite different 
values for a and b compared to previous works and this 
illustrates that the values of a and b are not independent 
given that the differences between the different works are 
not nearly as striking in Fig. [7] As a check of the use of 
extrapolated FM90 parameters to deredden the compar- 
ison stars (£ !2.1j) . we have verified that the linear fits at 
each wavelength do not change significantly even when 
sightlines with lower reddenings are excluded from the 
fit (up to A(V) = 1.5). 

We fit our a(x) and b(x) values using a similar func- 
tional form as given by CCM89. Given that we have 
seen in this work that the FM90 parameterization of 
the far-ultraviolet (x > 5.9 /^m _1 ) provides a good rep- 
resentation of extinction curves to i = 11 /im -1 , we 
have expanded the term used by CCM89 to represent 
3.3 /jm' 1 < x < 8 /im -1 to represent the entire range 
from 3.3 ^m" 1 < x < 11 ^m^ 1 . The term CCM89 
used for 8 jum -1 < x < 11 /im -1 was an extrapolation 
of the FM90 far-UV term. Thus, for the entire wave- 
length range considered in this paper (3.3 /im -1 < x < 
11 /im" 1 ) the functional form we use for our R(V) de- 
pendent relationship is 

a(x) = 1.896 - 0.372a; - 0.0108/[(.t - 4.57) 2 + 0.0422] 
+F a (x) (12) 

b(x) = -3.503 + 2.057.x + 0.718/[(z - 4.59) 2 + 0.0530] 
+F b (x) (13) 

where for (5.9 < x < 8) 

F a (x) = -0.110(2; - 5.9) 2 - 0.0099(a; - 5.9) 3 (14) 

F b (x) = 0.537(x - 5.9) 2 + 0.0530(a; - 5.9) 3 (15) 
and for (x < 5.9) 

F a (x) = F b (x) = 0. (16) 

We compare the average curves calculated from our 
newly derived R(V) dependent relationship with those 
given by CCM89, F99, and VCG04 in Fig. |1 The R(V) 
dependent relationship we derive is not noiseless and we 



show the uncertainties on the average curves to illustrate 
this point. The uncertainties on the curves resulting from 
uncertainties on a{x) and b(x) can be approximated by 
uncertainties on the value of R(V). Overall, these uncer- 
tainties are equivalent to approximately a 10% error on 
R(V). The uncertainties are significant and grow with 
decreasing wavelength. Overall, our new R(V) depen- 
dent relationship is roughly consistent within our un- 
certainties with past measurements in areas of common 
wavelength coverage. But, over most of the ultraviolet 
we do find a stronger dependence on R(V) than previous 
works. The statistically significance of this stronger de- 
pendence is difficult to assess given that previous works 
did not include an analysis of the uncertainties in their 
derivations of the a(x) and b(x). For VCG04, some of 
the differences can be traced to this study using un- 
weighted linear least square fits to determine a(x) and 
b(x). We checked this effect and found much shallower 
R(V) dependences with our sample if we perform un- 
weighted fits. Taking into account the uncertainties on 
both A{\)/A(V) and R(V) is clearly important. For 
CCM89 and F99, the largest differences are seen for 
x > 7 /im _1 where we find a stronger dependence on 
R(V) . This is not surprising as only with this new work 
has this wavelength range been systematically studied. 

3.4. Structure 

The existence of structure in the UV extinction curve 
can give strong clues to the carriers of different dust 
grain components. The main UV absorption feature 
is the broad 2175 A bump w hich has been i denti- 
ficd with carbonaceous grains (|Stecher fc Donnl H9651: 
iDraine fe Malhotral fl993l ). The other main structure 
is the far-UV rise and this feature is also identified 
(but with le ss con fidence) with carbonaceous grains 
|joblin et alJ Il992t iLi fc Drainel l2001h . Other than 
these two feature s, no other convincin g UV features 
have been found. IClavton et al . (2003a) has presented 
the most sensitive search for UV absorptions to date. 
They observed two heavily reddened sightlines and de- 
rived a 3<T upper limit of ~0.02A(V) on any features 
20 A or wider. Thi s resu lt has been strengthened by 
iFitzpatrick fc M assa (2007) who found a 3c upper limit 
of ^0.06A(V) on features 10 A or wider from the average 
residuals of 318 extinction curves. 

We have used our sample of 75 FUSE extinction curves 
to search of structure in the far-UV region. Fig. [10] 
presents the average extinction curve and residual curve 
for the entire sample at a spectral resolution of 250. 
There are two obvious features in the IUE spectral re- 
gion and these correspond to either ISM or wind line 
mismatches between the reddened and comparison stars. 
There are no obvious features in the FUSE spectral re- 
gion other than the scatter getting larger around x = 
9.25 /xm _1 and then getting much larger after x = 
10.75 /Ltm -1 . The scatter of the FUSE residuals puts a 3cr 
upper limit of ^0.12 A(V) on features with a resolution 
of 250 (~4 A width). For comparison to previous stud- 
ies, our 3<7 upper limit in the IUE range is ~0.04A(V) 
on features with a resolution of 250 (~8 Awidth). 

We can also search for features at significantly higher 
resolution given that the native resolution of our FUSE 
spectra are on order 10 4 . Fig. [IT] presents the average 
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Fig. 9. — The average extinction curves for R(V) values of 2.5, 3.1, and 5.0 are shown. The new average curves derived in this paper are 
shown as well as those from CCM89, F99, and VCG04. The uncertainty on our R(V) dependent relationship (i.e., a(x) and b(x)) is shown 
using the R(V) dependent fit that best describes the behavior of the uncertainties on the fits at each individual wavelength. In general, 
the fit uncertainties translate to an equivalent uncertainty of approximately 10% on R(V). 




Fig. 10. — The average extinction curve for the entire sample of 75 sightlines is plotted along with the average residual curve at a 
resolution of 250. The average residual curve was created by averaging the 75 individual observed minus FM90 fit curves. The deviations 
seen at 6.5 and 7.5 (im _1 are the result of ISM or wind line mismatches. 




Fig. 11. — The average extinction curve for the entire sample of 75 sightlines is plotted along with the average residual curve at a 
resolution of 10 4 . Overplotted on the average extinction curve is the H2 and H I absorption spectrum (dotted line) for HD 122879 
(A(V) = 1.41, R{V) = 3.17) to provide a guide to where strong H2 and H I absorptions occur. The average residual curve was created by 
averaging 75 individual observed minus FM90 fit curves. Overplotted on the residual plot as dashed vertical lines is the location of the H I 
absorption features. 
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Fig. 12. — The average extinction curve for 8.4 < A -1 < 9.0 (im _1 for the entire sample of 75 sightlines is plotted along with the 
average residual curve at a resolution of 10 4 . Overplotted on the average extinction curve is the spectrum of HD 122879 and the H2 and 
H I absorption spectrum for this sightline (A(V) = 1.41, R(V) = 3.17). The average residual curve was created by averaging 75 individual 
observed minus FM90 fit curves. The locations of features that produce residuals in the average extinction curve are labeled. 
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FUSE extinction curve and residual curve for the en- 
tire sample at a spectral resolution of 10 4 . There are a 
number of features in the residual curve. These features 
can be seen to correlate with strong absorptions in the 
example H2 and H I absorption spectrum. For exam- 
ple, two of the largest deviations correspond to Ly/3 and 
Ly7 at 9.76 and 10.29 /im" 1 , respectively. In addition, 
there is a strong residual at ~10.15 ^m -1 that corre- 
sponds to a particularly strong H2 absorption complex. 
In general, all of the strong residuals in the region be- 
yond 8.9 /1m -1 correspond to regions of strong H I or 
H2 absorption. There are residuals at x < 9 /im -1 that 
are not associated with H I or H 2 absorptions. Fig. [12] 
shows a blowup of this wavelength region including the 
spectrum of HD 122879 as a guide to the origin of the 
stronger residuals. The strong residuals all correspond 
to strong stell ar (C III, S IV, and P V), ISM (Fe II), or 
airglow lines (|Pellerin et al.ll2~002f ). Examining Fig. [Til 
the scatter of the FUSE full resolution residuals clearly 
becomes large where x ~ 9.6 /im -1 . Thus, the 3<r upper 
limits are ~0.15A(V) for x < 9.6 ^m" 1 and -0.68A(U) 
for x > 9.6 /im _1 on features with a resolution of 10 
(-0.1 A width). 

The average far-UV extinction curve can be compared 
with existing dust grain models to evaluate how well 
such grain models do at describing and/or predicting the 
far-UV extinction. The comparison is shown in Fig. [TO] 
where we show the average of all 75 sightlines (same as 
in Figs. [TO] & [UJ and the average of the 41 sightlines 
that have R(V) values within la of 3.1. The average 
R(V) for the full sample is 3.33 while the average for 
the R(V) = 3.1 sample of 41 sightlines is, not surpris- 
ingly, 3.13. In addition to the o bserved average curves, 
three different dust grain models (iWeingartner fc Drainel 
l200lHClavton et al.ll2003ll IZubko et al.ll2004D are shown 
in this figure. These models used different constraints 
and methods to determine the dust grain size distribu- 
tion and have known differences as a result. While all 
of the dust grain models shown are for an R(V) = 3.1 
(diffuse ISM), they are in closer agreement with the av- 
erage for our full sample than the R(V) — 3.1 average. 
This difference can be mainly traced to the differences 
seen between our work and prev ious determinations of 
the R (V) de pendent relat i onship . iWeingartner fc Drainel 
(|200lh and IZubko et all (l200l used the R( V) = 3.1 
curve from F99 while IClavton et all ((2003b) used the 
R(V) = 3.1 curve from CCM89. As both the CCM89 
and F99 R(V) = 3.1 curves are below our R(V) = 3.1 
determined from the R(V) dependent relationship (see 
Fig- ©j it is n °t surprising the dust grain models are also 
below our curve. The differences between the CCM89, 
F99, and this work are due to different samples of extinc- 
tion curves used as well as dif ferent fitting appro aches. 
The lClavton"et~aII (|2003bft and lZubko et~aLl (|200l mod- 
els do a much better j ob in reproducing the far-UV ex- 
tinction lev els than the |Weingartner fc Drainel (]2001l ) (as 
updated bv lDraing|2003f ) model. Overall, the far-UV rise 
seems to consistent with its potential origi n as a second 
resonance feature center ed around 14 /im _1 (poblin et al.l 
[1991 iLi fc DTainl[200l . 

4. CONCLUSIONS 

We have presented FUSE+IUE extinction curve for 75 
sightlines that have A(V) values from 0.6 to 3.2 mag. 



The R(V) values of this sample range from 2.5 to 5.5 with 
a strong clustering around 3. These extinction curves 
were created using the standard pair method. Given the 
FUSE sensitivities, this required generating a new set of 
unreddened comparison stars. This set of 75 FUSE ex- 
tinction curves represents a large increase in the number 
of extinction curves well measured in the far- ultraviolet. 

Using this sample, we investigated the nature of the 
far-UV extinction with three methods. The first was to 
fit all the curves with the FM90 parameterization. We 
found that the FM90 parameters determined with just 
the IUE portion of the extinction curve to be broadly 
consistent with those that were determined using the 
FUSE+IUE extinction curves. There were significant 

changes in the individual values of the C^ V ' and 7 pa- 
rameters. In addition, the average parameter de- 
creased by —8% indicating that the strength of the far- 
UV rise is overestimated when using only IUE data. All 
the extinction curves had positive far-UV rises when fit 
with the FUSE+IUE data giving no indication that the 
far-UV rise was turning over at the shortest wavelengths. 
We tested the new FM07 parameterization (with 1 addi- 
tional parameter and a simplified far-UV rise term) using 
the combined FUSE+IUE curves and found only a weak 
preference for the FM07 parameterization. We choose 
to use the FM90 parameterization as it is simpler and 
provides for direct comparisons to previous work. 

The second investigation of the far-UV extinction con- 
centrated on deriving the R(V) dependent relationship. 
We derived the linear dependence of A(X)/A(V) on 
i?(U)~ 4 at each FUSE and IUE wavelength including 
the observational uncertainties as part of the fitting. We 
found a stronger dependence of A(X)/A(V) on i?(V) -1 
than previously indicated in the far-UV. It is not sur- 
prising as the previous work in the far-UV was based 
on only a very limited number of sightlines. We in- 
clude the uncertainties in the fit in the R(V) depen- 
dent parameterization and find that ou r derivation is 
broadly consistent with previous work ( Cardcll i et alJ 
H989tlFitzpatrick fc Mass jl99fl IValencic et al.ll2004f ) in 
the regions of overlap (mainly the IUE range). We 
present our R(V) dependent relationship for the entire 
UV (3.3 iixnT 1 < A -1 < 11 /im -1 ) using the same func- 
tional form as used by CCM89. 

Finally, we searched for discrete absorption features in 
the far-UV. We looked for features at a resolution of 250 
by averaging all 75 extinction curves and also averaging 
the residuals of each curve from its corresponding FM90 
fit. We find a 3tr upper limit of ~Q.12A(V) on features 
with a resolution of 250 (—4 A width). Utilizing the 
high spectral resolution nature of the FUSE observations, 
we also searched for features at the native resolution of 
10 4 . We found 3cr upper limits of -0.15A(U) for A -1 < 
9.6 jam" 1 and ~0.68A(U) for A" 1 > 9.6 ^itT 1 on features 
with a resolution of 10 4 (—0.1 A width). 

We thank Derek Massa for useful conversations on ex- 
tinction curve uncertainties and Karl Misselt for discus- 
sion of dust grain models. We appreciated the comments 
of the anonymous referee that resulted in a stronger 
paper. This study was supported by NASA ADP 
grant NAG5-13033. Based on observations made with 
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FlG. 13. — The average extincti on curve for the 41 curves in t h e sample that have RC V) va lues within la of 3. 1 is shown at a resolution of 
250. The R(V ) = 3.1 models fromlWeingartner fc Drainel Hwmb , [Clayton et al.l J2003bT ), and lZubko et ail <200lD are plotted for comparison. 
Note that the Clayton et al. (2003b) model only extends to 10 (an -1 . 

the NASA-CNES-CSA Far Ultraviolet Spectroscopic Ex- kins University under NASA contract NAS5-32985. 
plorer. FUSE is operated for NASA by the Johns Hop- 
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